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I. INTRODUCTION 
[iieE@e Js am Lnezeaesed emphasis in the ability ¢o fcza- 
cast the behavior and evoluticn cf aerosols. The sg 
and aksorpticn of light are influanced by aerosol 
tions in the atmospheric boundary layer ane a S 
perfcrmance of optically guided weapon systems. The Air 
Force is particularly interested in how aercsol extinction 
affects the use cf precision guided munitions (PGM) 
kot=reli et al, 1979/ The Department of Defense (LOD) has 
PGM's that cperate at differing wavelengths which range from 
the visirkle to tke microwave regions. PGM's have a greater 
meeeaaty tc hit a target than conventional munitions; 
however, the eenrreliliaag factor 1S the abilivy of -the 
guidance system to “see" the target. This ability is deéepen- 
dent cn the wavelength for which the sensors are designed 
and the properties in the intervening atmosphere. The 
degrading froperties in the atmosphere are principally 
tolecular absorption and aerosol scattering. The wavelengths 
for tke quiéance systems are designed such that the mclec- 
ular abscrption is minimized: therefore, scattering Ey aero- 
sols becomes the main concern once a suitable 


absorpticn-free window has Eeen selected. 


The model under ccnsideraticn includes the behavicr of 
fMarine aerosol compcnents as well as previous continental 
compenents. Estimating the influences of marine aerosols on 
electro-cptical (£0) systems has been studied [Barnhart and 
Srceete, 1970/ » Particle sizes of interest are those asso- 
ciated with locally generated sea salts because cf their 
effects on IR as well as visible wavelengths. The size 
distributions of sea salt particles show @ variance of 


several crders of magnitude. 





f[eeematodicy to Z¢recest the behavior of aerosol ex*tizc- 
tion frem synoptic scale patterns would help in the dec 
which tyre cf weapon system to employ. Because scme systems 
pmemeauncred Ercm tke air, it is important that eu 


Baegetaile include the vertical distributions of aerosol rerti- 


O 
cles. Weodels ss weckiSt “or  eStimacing “vertical  extincticr 
beonales, Eut they have not been sufficiently verified. Teo 
do this, profiles cf actual aerosol data amu 


and ccmpared witk the model fcrecasts. 


Mcdels in current use are based on paramsterizaeticns of 
the effects of relative humidity and wind speed on the equi- 
ieeeaiut aercsol distributions jaells, et aL. 1977/ s- RECENT 
evaluations have shown that these models are limited te mean 
distributicrs (i.e., the average aerosol concentration ata 
Given wind speed and humidity) (Padsal1 etal, 198 2a / ; 
Models are limited because scmée processes in the atmnespheric 
mixed layer which affect aerosol concentrations are not 


considered, namely entrainment and inversion height changes. 


he purpese of this study is to present and evaluate 4 
model which includes the meteorological processes that 
adequately describe the whele marine atmospheric tLoundary 
layer (MABL). An inversion represents a cap to the vertical 
transport of surface generated aerosols, 2nd. -uS- noe 
accourted fcr in previous models. The top of the bcundary 
layer is capped by the marine inversion, where entrainment 
cf overlying air takes place. Because €ntrainment mixes 
clear (ncn-marins) air inte the marine layer, this precess 
1s aS important as surface layer fluxes in determining equi- 
Peoria un ccncentration. This entrainment process is included 


in the mcdel. 


ite 
from ar experiment entitled Marine Aercsol Generati 
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Et (MAGAT). The experiment was conducted in tthe 
feernoxy Of Monterey Eay, Ca, during the perioc cf 28 April 
Semo fay 1980. The purpese of the experiment was to examine 
Grewce@ratibility of cptical and micrometeorologicai crcopfa- 
Geelon =heory, and tc extend dynamic models of the evelving 
MABL ele include aerosol and turbulence bao es 
[Fair all, 1980 and Fairall et a1,1980/ a Pawo ble SE cams 7s cic 
R/V ACANIA, and an aixrcraft were used. 


Meeents Study, twe 24-hour periods were chosen for 2valu- 
meme. 3 May, beginning at 1200 PDT and 6 May, beginning at 
feo ELT. These pericds were chosen beacause the boundary 
layer was undisturbed (no fronts closer than 100 nm) fer 24 
hours pricr to the starting times. Additionally, the times 
were sélected dte te the proximity of the aircraft and 
surface shir during the experiment times. The approach was 
to describe the synoptic conditions from 24 hours in advance 
cf the mcdel forecast through the end of each forecast, and 


compare the model output with the actual MAGAT findings. 


11 





If. BACKGROUND 


ae oe ee ee ee ee 


A. DESCRIFTION OF EXISTING MCDELS 


Currert sodels for estimating aerosol equilibrium distri- 
tutions use meteorolcgical inputs of (10-m height) wind 


speed ard relative humidity. These two quan 


mw 
23 in 


ie 
censidered tecause of their rele in generation end transport 
(wind) and growth (humidity) of aerosols. The Navy' 
woe -Munn=Katz (WMK) fwelis = Oeil, 1977] is an example of 


mas CCNCept. The performance of this typa of model has 


= 


meen Studied with data obtained in the northern Atiantic and 
eastern Facific Ccean areas. The model output compared with 
these data is shcewn in figure z.1, depicting a height depen- 
dence cf tctal aeroscl volume from a sampie set of eastern 
Facific data. These profiles correspond to (1) the ckserved 
sea salt vclume, V3 (2) the observed sea salt veclume 
adjusted tc 80% relative humidity, VW; and (3), the WMK 
predicted volume adjusted to 80% relative humidity 
(circles). it is clear frem this figure that within the 
mixed layer the observed decrease of aerosol volume with 
height is less than the model predicts. The surface geéener- 
ated aercscls appear to be well mixed below the inversion 


when normalized to 80% relative humidity. 


The assertion that existing models can predict only a 
mean value appears in the results obtained in the North 
Atlantic, from the JASIN Experiment. Figure 2.2 compares a 
Single radius size cf five micrens and the corresponding 
model prediction. The values and trends in the predicted and 
mean results are in reasonable agreement. However, the 
Standard deviaticn is three times the mean; if one assumes a 


normal distribution, cnly 67% cf the observed aerosol 
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Figure 2.1 Height dependence of ASR CSosy volume and relative 
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Fagure 252 Aerosol volume spectrum at r=5 microns 
as a function of wind speed, 
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Pesce iteuclons will ke withir a factcr 
average. This comparison emphasizes the o 
how accurately a mocel predicts an évera 
at a given wind speed and relative humidi 

c 


merece RMS veriations cannot ke elimineai 


J 


€mang mcre reteorolcgical parameters. 


EF. THE NES BOUNDARY ILIAYER SODEL 


{he distinguishing feature of this modél is the charac- 
memezacicn cf the MAEL, which is convectively nixed up te a 
Fergnt h, and capped by an inversion. The atmospheric 
profile is depicted in figure 2.3a, representing a cloud- 
free mixed layer where water vapor mixing ratio (g), and 
Virtual fotential temperature (e,) ace "well-mixed", ie, 
independent of heigktt in the mixed layer. The assumed 
vertical aerosol profile is shown in figure 2.3b. The model 
produces 2 24-hour time evclution of an aerosol Seectsun, 
requiring a predicticn of the fcllowing at each time step: 

1) surface producticn rate of marine aerosols 

2) entrainment rate at the top of the mixed layer 

3) mixed layer depth. 
Tetails cf how these parameters are input into the model are 
discussed in chapter IV. The mathematical relationship 
retween the time rate of change of the aerosol volume spec- 
ee,  avV/dr, and these three parameters is shown in the 
follcwing equaticn 


dv/dt =(<W'V'>o - (WetWkm) v7 A 


The mcdel Seaqwere. the €<VOlUtICH sOf aSsrosol at. ‘five radii 


meee ae oso, 10.05 \and 15.0 microns) of both the ccnti- 


nental and the marine (sea salt) components in the mixed 
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mir, Lhe ectual tine evoluticn cf each radius is cutrut in 


the fcrm of dv/fdr, with the units of um? /cm3. 


Dae modei assumes ro appreciable concentraticn of wearin: 


(Pp 


aeroscis akcve tke mixed layer. The entrainment cf air i: 


by 
ft 
O 


3 
mv 
rf 
t-4- 
+3 
(D 


the mixed layer will not caus¢ an element of thes 


¢ 
particles tc "escape" into the free troposphere. ane §ce 


ci 
sv 
J 


ct 


h 
h 


Marine ccmpcnent. The ccean also acts as asink fer “he 


acts as a scurce for marine aerosols, primarily through 


iD 


cr 
mM 


t 


Bememecticn cil White caps. This is the only input for 


h 


continental aerosols, which are generally cfa smaller 
a 


S 
radius. The entrainment precess at t top of the mixed 
REO = hes AEE. Dass 


process could be as important as the surface layer fiuxes in 


dayer mixes clear (ncen-marins) ales a 


deteéergining an aeroscl concentration. 


Cw. SMCDEL INPUTS AND AEROSOL INITIALIZATION 


The model is designed so that with the exception of 
initializing the aercsols, all of the calculated parameters 
are kasec cn inputs from surface based cbservaticns, and 
representative soundings. The input consists of: 

1) Non-meteorological inputs of latitude, julian day, and 
the leccal start time. These are used for a diurnal radia- 
tive heating/ccoling package. Diabatic warming has an 
obvious impact on the life and strength of the rixed 
Mayer. Additionally, sea-surface temperature (SST) is 
input at the start time. Up to ten forecast SST values 


Can ke infgut during the 24-hcur period as well. 


2) Surface wind speed and direction at the start time and, 


as with SST, up to ten additional forecast values. 


3) Frem figure 2.3, temperature and relative humidity are 


parameterized by virtual potential temperature and water 
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Maer Wiking ratio, respsactively. Mixing ratios are input 
above the mixed layer, as well as the lapse rate 
oa p! 


mee Ne tializa, aon or (thea. vsix- > radii cencentra 
Mmpee Lecires is reguired. To do this it is necéssa 9 
understand the actual distributions of the continental and 
Marine ccmpenents in the mixed layer. For «he purposes of 
ip 


€valuating the model, assumptions must be madé¢ cn ?¢ 


oe 


h 


Mes raputicnrs of the acquired WAGAT data. The aircraft 
Mimecrcuments could nct distinguish the chemical make-up of 
miemandividual particulates; therefore, it is not knewn, for 
aeotven radii concentration in the mixed layer, hew much is 
continental and how much is marin2. The concentraticns of 
€ach size were caiculated in the mixed layer, and then abcve 
the mixed layer. Based on the assumption that there are no 
marine ccencentrations above the mixed layer, we conclude 
that the concentraticns calculated from abcve are all centi- 
rental, and therefore linerally subtracted from the concen- 
mae 2 Cons in the PExed Selayver, leaving only marine 


G@oncertraticns. 


S 


f-4¢ 


A further adjustment of the input a¢grosol values 
required, tased on changes in relative humidity during the 
data collection time. We have stated that the growth of 
aeroscls is a functicn of relative humidity. Therefore, a 
distributicn of aeroscls gathered at 90% humidity cannct be 
directly cerpared to another distribution gathered at 80 % 
humidity; 2a reference humidity is required. Consider the 


aerosc] welume spectrum 


vir) = 4/3 er3 n(r) (2.2) 


ie 





where n(zr) is the nugker density spectrum. V(r) ie ae Cer 
as the vclume of aercsol particles per cm? at a referen 
Eeeuacicn ratio of S=80%. A humidity growth f£ 
fFairall, eval, 198 2a/ is defined as 


G(S) =.81 exp(.066S/(1,58-S)) (2.3) 


There are further considerations concerning advecticn. 
mleaciy ell aerosols do not originate lcecally; both above 
and below the inversion eeroscls are advected inte a local 
regicn. In terms of aerosol density, entrainment acts as an 
aeroscl flux out of the boundary layer because the ccncen- 
traticns abcve and below the layer area different. In the 
fhemel, the entrainment acts on the “jump” across the inver- 
S0n. For the purpcse of €valuating the the evoluticn of a 
local ccncentration cf aeroscls, the medel will neglect 
horizcntal advection and further assume a negligible local 
producticn cf the centinental ccmponent. The mcedel takes 
into ccnsideraticn a Stokes gravitational fallout term W, 
fan, 1679] - The fallcut rates above and below the inversion 
are different because of the changes in aerosol spectra, 
caused by the humidity growth factor. The Stokes velccity 
frairall, G2 ~ aly 16 &2a/ is calculated from the fcllicwin 
eguaticn where p, is the density of the droplet and E is the 


kinematic viscosity cf the air. 


2.4 
W, = 2g(p, -p) r2 G24S) / (9Ep) eee 


uo 





CT. EQUILIBRIUM ABEROSCL SODEL 


Meee a Means Of  Comparis¢cn cf asrosol initializetion 
schemes, the model cffers the option of initialization with 
equilikrium values. These initial valu 


collected from the JASIN Experiment 


+- 


re) 


S are bas2da on data 


tz] 


rall, e+ al, 1982/ . 


ad 
=> dle ny -equii st rroun 
2 


A large aercsol data base was collect 
@emarticrs, with respect to wind snvneed, Giese Gg. Zee he i 
periods and then nermalized to a reference relative 
numa dat y. Aerosol spectra were then grouped into six 
gifferent wind sreed rangés. Fer the equilibrium initializa- 
tion scheme in the NES model, a pre-assigned value cf «he 
five radii are used as initial aerosol values, based on the 
wind speed at the mocéel start time. An example of the equi- 
libriur aercesol spectra from the JASIN Experiment is shcewn 
mec agure 2.4, Note tkat the reference relative humidity for 
this experiment was &7 percent. The graph is referenced here 
See tc chow an equidibrium type distribution; the initial 
equilibrium and MAGAT aerosol values will be shown in 
Gmapcer IV. Once initialized with either equilibrium data or 


calculated data, the model will run in the same manner. 
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Figure 2.4 Aercsol spectra as a function of size 
Ox select ed wind speeds. 
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This ¢quilibrium scheme tliustrates the fundamenczal 
difference between the NPS model, which is a Gynamic medel, 
gid a Sc¢ady-state tyre model. Ak representative me 
mags eguilibrium approach 2S the WMK model, mencicnee in ch 
Beqeannane Of this cherter. This type of output describes the 
Werane aeercsol distributzons as a furction of wind speed 
(surface generation), relative humidity (growth factor), and 
elevaticr {vertical variation of aerosols with height 
assuming a steady state vertical transport pire 
mathematical representation of this model [Faizall, e. al, 
1982a/ is 


y 
n(r) = (r/a)+1.62 (Cy+CoV )F*eXp (-2/hgF-8.8(r/a) (2) 


where 
moet hee Part mele: saze =n Miecsons , 
u= tke wind speed, 
v= 0.5 for us 4 Ps, 
VerUe—-35.5)- 105 U0 4 0/7s, 
f= 1+(v/60) 3, | 
Y= 0.364-0.00293v'29 , 
Z= keight above sea surface, mn, 
heeewscale height,n (800 me for 2<1000 no), 
a= 0.8lexp (0.066s/(1.058-s)) , 
S= Hy/100, (H= relative humidity in percent). 


The cther constants are given, based on the value cf v. 


Voss e. Co 6 
was 7 350 1000 West 
v>7 0 6900 aco 


The equaticn described produces a number density spectrum, 
n(r). This relates tc the volume spectrum V(r), as shewn in 


Colac dOn Ze Zs 





III. DATA ACQUISITION AND SYNOPTIC CONDITIONS 


Aw CATA ACCUISITION 


The cata were obtained frem a ship and aircraft exreri- 
ent, Marine Asrosol Generatica and Transpor= (MAGAT), held 
maemeeco APril to 9 Mey 1980, in a region 30 to 50 miles cif 
fmameeeccasSt C{f Monterey Bay, Ca. Data gathering was dere with 
instruments mounted cn both an aircraft and é&@ surface ship, 
the RyV ACANIA. Aeroscl data were obtained from the Airkorne 
Research Associates turbo-charged Bellanca, using a Farticle 
Measurinec System (PMS) Axial Symmetric Scattering Aerosol 
Probe (ASSAP) rarticle counter. All measured data were 
sampled every 2.5 secends, with @ two-scan average cf every 
muverceccnds. The scans were ccllected in “ladder" profiles, 
during which the aircraft made measurements at a constéent 
mer tuace £cr two minutes, climbed to a new altitnde, and 
repeated a rew measurement run. The instrument utilized 60 
Size channels from 2.8 to 14.0 micron radius. In most cases 
the ladders extended from near the sea surface (3 m) up up 
through the well-mixed boundary layer, to a few steps above 
the inversicn. The elevation was generally up to 5 kilome- 
ters. A typical ladder profile contains 10 to 14 steps. The 
step heights were randomly chosen, but an attempt was made 
to keep each step height consistent between ladders. Air and 
dew-rcint temperatures were also measured and used te calcu- 
late relative humidity for the correction factor menticned 
in chapter two. The aircraft also flew ascending spirals (in 
the vicinity of the ladders) during which other metecrolo- 
gical parameters were collected. This data yielded vertical 


soundings similar to those pfrevided by radiosondes. 


eee 





There were two aerosol instruments fitted on tne R/V 
ACANTA. One was the PMS model CSAS (classical scattering) 
and the cther was model ASAS (active scattering), ccntrclied 
ky a EMS data acquisition system (DAS-32) With) a COT uces 
ace rface. The shirkoard system measured aerosols in 90 


@@mEnrerent size channels from 0.9 to 14.0 micron radius. 


It was nected, during the flybys over the R/V ACANTIA, that 
the aircraft measurerents did not agree with that of the 
Ship. The size distributions from the aircraft measurements 
were consistently smaller than the ship measurements, for 
madeg=eater than 1.0 micrcn. In addition, the differences 
miereaced with radius. Since the shipboard aerosol system is 
newer and hed a wider size range, and better sensitivity, 
aircraft aerosol data were corrected to agree with the shi 
aerosol Gata /Fairall,1980 and Fairall et al,1980/) . 


Prefiles of virtual potential temperature and mixing 
ratic were cbtained from three different sources. The reason 
for usinc these parénmeters instead of temperature and dew 
FCoint, is that mixed layer inversions are more easily idén- 
tified with the former variatles. The sources were the 
spiral flights from the aircraft, the radiosonde launckes 
from NES and from the R/V ACANIA. 


E. SYNOFTIC DATA 


Surface and 500 ab syncptic charts and the GOES WEST 
Satellite images were used to evaluate the synoptic ccndi- 
jm. On S « Charts are from the NOAA weekly series cf daily 
weather maps. In addition to data collected from the R/V 
ACANIA and the aircraft, local weather data were also avail- 


able fier the U.S. Army's Fritzsche Field weather facility, 
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Figure 3.3 Same as Pigure 3,2 except 1245 PDT, 3 May. 
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Figure 3.4 Same as Figure 3,2 except 1645 PDT, 4 May. 
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Figure 3.5 Same as Figure 3.1 except 5, 6, and 7 May, 1980. 
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Figure 3.6 Same as figure 3,2 except 1745 PDT, 5 May. 
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Figure 3.7 Same as figure 3.2 except 1645 PDT, 6 May. 
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Figure 3.8 Same as Figure 3.2 except 1745 PDT, 7 May. 


nt 











bere Crd, Ca. Satesiiee EnCtecdraphs and NOAA mape ter the 
two different time periods (2-4 May and 6-7 May) areé she 
in figures 3.1 through 3.8. Each three- day sequence dericts 
the ZU tcur period t¢ing censidered; 24 hours b¢fere ane 
Seiten. Atmcspheric soundings for the two model runs are 


Efewn in chepter IV, figures 4.1 and 4.2. 


C. SYNOEFTIC CONDITICNS 


t+ 


Several weak frontal systems passed through the area 


aang the experiment. Showers occurred during the fiztst an 


Vj 
wy) 


$b Oy 


tase Gays, associated with the fronts. Lew cloudiness an 
mecmecccured during the morning from 29 April to 5 May, with 
ESgeLeEturnikg again cn the last day. 


At the beginning cf the period, the area was dominated by 
a slcwly eastward micration of a cut-off low at 500 mk. By 
eemey mMCrIning on 2 May, the area was under the influence of 
aweak ridge. On 3 May the area was under divergent flcw at 
the upper level. An upper level low had formed off the Baja 
Ca. ccast cn 4 May, leaving the area under the influence of 
acol. Cn 5 May the area was between an upper-level trough 
and ridge, and by 6 Fay the area was on the back side of the 
trough. Eecause cf tke deepening of th2 trough, the area was 
esoemiet CX the back side of the trough on 7 May. A new upper 
level trcugh formed énd appreached the area on the final two 


days cf the experimerc. 


Surface winds were relatively light, 0 to 10 kts, at the 
beginning of the pericd, and increased toward the end of the 


Pemted tec 16 kts, gusting te 22 kts. 


An inmpertant feature of these interpretations is the 


nature of the mixed-layer, cften topped by an inversicn, 


c 
With reguard to stability, and therefore mixing intensities. 
It is assumed that tte mixing becomes greater as conditicn 


tecome mere unstable. 





Mae itemavd= cs €h cle e€xXpetament, the mixed Jayer exhibited 
maeeee te Sflightiy unstable conditions, until abcut 1800 
Mecal tine (PDT) 28 April, when conditions became mere 
moweradal, The neutral conditicn remains until 1 May when 
conditicrs cnce again become stabilises. A weak frontél pa 
memore 05S00 PDT on 2S April does not appear tc affect the 
mixed layer profile. The layer remains staktle until a 
frontal passage on 2 May when conditions beccme neutral and 
Femairs sc tntil 5 May. On the morning of 5 May, conditicns 
aze Slightly stable, PG Eetunrned LO Tneltral “on, 6 May, 
zi 


despite a frontal passage at 1300 PDT May, and 


OA 
remaining neutral to the end of the experiment on 9 May. 


C. DATA SELECTICN 


Data from ths MAGAT Axperiment were selected for model 
verificaticn on the kasis of the relative pcsiticns of the 
meo@mscata Gathering platfcrms and the general synoptic ccndi- 
tions. The locations of the surface ship and the aircraft 
did not always coincide during the experiment. Therefcre 
ecundirgse and aerosol data were considered suitable if the 
two were within 20 nautical miles. Hourly wind srpeed, 
direction, and SST were gathered from the ACANIA and were 
used alcng with th¢ aircraft soundings. PLOts CL. “wind 
speed, air (solid line) and sea-surface (dashed line) temp- 
€rature, and relative humidity for the two 24 hour periods 
Beewociown 2n figures 3.9 and 3.10. Another consideraticn of 
the aircraft data was the relative locations of the aerosol 
ladders and the soundings. Although both gathered during 
vertical flight profiles, they were not done at the same 
time. Tke synoptic ccnditions during the experiment include 
pericds cf frontal weather which was alse avoided 1 
analysis. Two time blocks 
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Figure 3.9 Plots of wind speed, air and sea surface 
temperature, and RH beginning 1 DT 3 May. 
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Figure 3.10 Same as figure 3,9, except beginning 1800, 6 May. 
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Figure 3.11 Flight eae of taitt on themmorning cf 
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ay, 1980. 
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Figure 3.12 Same as figure 3,11 except afternoon of 3 May, 1980. 
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Figure 3.13 
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Flight path of aircraft on 4 May 1980. 


So 





FLT #12 oa 





6 MAY 80 
37 ON 

SANTA 
CRUZ 

mee ge 

jesse 
18So 
1@sy 
IfSe 
SE oN 
0 10 20N.Mi 122%, 


Figure 3.14 Same as figure 3.13 except 6 May, 1980. 


52 





FLT #13 





f MAY 80 
S7 Hy 
SANTA 
Mee S CRUZ 









CS 






}Oo| 
Seas 





T a more en 
\ SA 

RY i 1§ 
| 


36°, 


Oo 10 20N.Mi. 
ee ee ee 


122) OW 


-” 


Figure 3.15 Same as figure 3,13 except morning of 7 May 1980. 
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Figure 3.16 Same as figure 3.13 except afternoon of 7 May 1980. 
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Figure 3.17 Cruise track of R/V ACANIA 3 May 1980. 
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Figure 3.18 Sage as figure 3.16 except 4 May 1980. 
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Figure 3.19 Same as figure 3.16 except 6 and 7? Say 1980. 
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Meme Chocen; 3-4 May and 6-7 May. Figures 3.11 through 3.16 
Mmgicats the routing of the aircraft flights for «hs two 
meme Ftlocks, and figures 3.17 through 3.19 depict the ccurss 


cf the RyV ACANIA during the same cimes. 


MemGsEneric data from the aircraft soundings were used 
instead cf the radicscndes frem the R/V ACANTA. This was 
done fcr twe reasons: first because the aircraft scundings 
containec mcre lévels, and second they were in closer rcrcx- 
Mery tc the aircraft aeroscl ladders. <tmospheric and 
wemerscl! €ata £or the first run wes taken from flight 6, cn 
the mcerning of 3 May. The locations and times are shown in 
figure 3.9. The aircraft was approximately 45 miles frem the 
coast when a soundino was taken at 1143 PDT, and an aerosol 
dJadder (111) was taken at 1200 PDT. Data for the two verifi- 
cations of atmospheric parameters and aerosol values were 
chosen tased on their proximity to the initial data. A 
sounding and aeroscl ladder cn the afternoon of 3 May 
(sounding time of 1714 PDT, and aerosol ladder 113) approxi- 
mately eight hours frem the initial time was chosen for the 
fist VETificaticn, and is depicted in the route cf flight 
feeeagure 3.10. Data for the second verification was chosen 
Deems tlight 8&8, GUGUure terials “cOuUn@iIng tate io lePC?.. vane 
aeroscl ledder 116. fata for the two verifications are beth 


wmaeenin 2C giles of tke initial data. 


Data for the seccrd model run is shown in figures 3-12 
through 3-14. The medel was initialized with data frem the 
1732 EDT scunding and ladder 122, approximately 25 miles cff 
the Menterey coast. The first verification data was taken 
from flight 12 the next morning approximately 17 hours 
later; sounding time 1743 PDT and ladder 123. The seccnd 
yverificaticn data was taken from flight 13, approximately 24 
hours frcem the initial data; scunding time 1836 PDT and 
ladder 124. Data ccllected for the two verificaticns were 


Within S titles of the location of the initial data. 
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Maqdvtaenal wind speed, a and SS? colle€a@aad f-=o0n 
n 


cf the tso model runs. AS mentioned in the model dsscrip- 
fen, ub <C 10 forecast values for each can be input into 
the mecdel. However, due to the continuous movement of the 
R/V ACANIA during tke MAGAT experiment, where SST values 
were changing because of strong coastal gradients, a single 
value of SST was used. This value was chosen when the B/V 
ACANIA was closest tc the areas where the aircraft data was 
meme ©rcr €ech cr the two initialization sites. The wind 
speed and direction cata freon the R/V ACANIA was considered 
representative of the data collection sites, and therefcre 
a 


CL 
1s. 


(D 


memo into the mcdel runs at arfproximately 3 hour interv 
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Aw DIGITIZED ATHOSPEERTC SOUNDINGS 


The atmcspheric irputs for the model (chap 2) azé: 
1. The mixed layer equivalent potential temperature (C) 
ang specific humidity (gm/kg). 
Zee the "jump" disccrtinuitiés in each of the above values. 
The jump meaning the difference etween the mixed layer 
value and the value at the tcp cf the inversion. 
7. The lapse rate for each parameter value akeve the 
Bmiversicon. 
Beeeuhe déerth of the well mixed-layer 


To simplify the frocess of calculating these values from ths 


A 


soundings, a digitizing scheme was designed by “Gh 
Environmental Physics Group at the Naval Postgraduate 
eemeci. ‘This scheme transforms the sounding into the struc- 
ture that descrikes the necessary model inputs. Examples arse 
shown in figures 4.1 and 4.2. Note that the "jumps" are 


@emerected as occuring in an infiritely thin layer. 


Ew ALDDITICNAL MODEL INPUT AND ADJUSTMENTS TO THE DIGITIZED 
SCUNDBINGS 


The reascn for making adjustments to the atmospheric part 
cf the mcdel is to previde a better basis for evaluatina the 
rehavior of the aeroscl prediction. This is optimized when 
the mcdel is preducing the best possible forecast cf the 
atmospheric parameters. The parameters are forecast very 
noe in Foth model runs, with these adjustments made. The 
model is accurately éescriking the height of the inversion 
through adjustments in the subsidence rate, and the lifting 


condensation level (LCL) is accurately generating clouds 
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through adjustments tc the temperature and specific humidity 
which are within the limits of instrumentation ertrer. ‘ea 
cloud verificaticns were made with satellite data and ctser- 


Mees=cnhs trem the R/V ACANIA lca. 


ME initial “#nput for the model is the subsidesce rate. 
Using a first guess cf -.0C5 m/sec, the model generated a 
Flot cf the inversion height for z=he 24-hour period. FEased 
on the verification data, the subsidence rate was further 
adjusted tc bring both the predicted and cbsarved inversion 
keight values together. Then a small adjustment was made to 
the mixed layer values to match the stability and cloud/ 
Cloud free patterns with what was actually observed during 
m@mempezticd. In the first model run (3 May), ‘the petential 
temperattre was increased 0.5 degrees and the specific 
humidity was decreased 1 gm/kg. The resultant suksidence 
rate necessary tc match the predicted and cbserved inversion 
hFeight was -.0042 m/sec. These adjustments yield an atmos- 
pheric profile similar to what was actually observed. The 
alr mass became slowly saturated enough to generate clouds 
in the upper part of the mixed layer at about 0300 PDI on 4& 
May. Ficts of the inversion and LCL behavior, and the temp- 
erature and humidity rfrofiles are shown in figure 4.3. In 
emo rlot fcr LCL and inversion height, ene -clettes "HH" 
depicts the verification times and values for the inversion 
height, andthe letter "L" similarly indicates values for 
meemncl verifications. The letters "fT" and "QO" in the temp- 
erature and specific humidity graphs indicate verification 


values fcr e€ach respective value. 


The atucspheric flets for the second nodel run (6 May) 
are seen in figure 4.4. For this run the specific humidity 
was reduced from 6.4E& to 6.3 gmyvkg, the potential temrera- 
ture was unchanged, and a subsidence rate of -0.01 m/sec was 


muput for the run. 
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Figure 4.1 Digitized soundi 
May 1980, 11 
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FOr first Bodel run, 
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Figure 4.3 Plots of inversion and LCL heights, temperature, 
and humidity profiles, first run. 
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Figure 4.4 Same as Figure 4.3, except for second run. 
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C. AEROSCL INITIALIZATION PROCEDURES 


The aerosol data were gathered in ladders, usually up to 


Fave kiicmeters, and were of actual concentrations cf parti- 
cles for a given increment cf radius, n(r). This number 
contained bcth marine and continental components, ccllected 


a 
at the anbkient relative humidity. 


The initializaticn scheme was perfcrmed in three main 
Beeps. The first ster was to identify the marine and centi- 
hnental ccmpcnents. By using the digitized soundings for each 
model run, & detergination was made of two height values: 
cne representing the well mixed-layer and the other reore- 
senting the free tropesphere, abcve the inversion. Both the 
tBixed layer and free troposphere aerosol counts for ¢ac 
the five radii were calculated. The second step was to aprfly 
the humidity growth factor +o the values and calculate 
av/dr. The relationship between n(r), which was measured, 
Gmeeav/dr 15 given by equation 2-2. The final step was to 
Subtract, fcr each réecius, the free troposphere dv/dzr valus¢ 
from the mixed layer dv/dr value. With the assumpticn of a 
constant prefile of continental aerosols through tk;oth the 
mixed layer and into the free treposphere, this subtraction 
leaves in cnly tha mixed layer marine components, normalized 
to 80 percent relative humidity. The values of dvyfdr used 
to initialize the mecdel, both equilibrium and MAGAT, aré 
Shown in tables 4.1 and 4.2. 


This three step precess was dene for both model runs, and 
for tsc verification times within each run. As described in 
Chapter II, the verificaticn ladders were chosen to Le as 
close as pcssible tc the position of the initial a¢rosol 
data. In koth cases, the verification data were within 15 


tiles cf the initial data. The plots of the aerosol output, 


a 





meri 6G the MAGAT and equilibrium initializations, aS 
Sepected aS logi0 (dvyvdr) against the 24 hour time scerzicd, 
Mmempiadlc Hour time sters. In beth (initializatior) cases, the 
Peeevalues represent the MAGAT verification values. In beth 
feet cUnS, aeroscl output for the MAGAT initiaibzetion 
schere is represented by a continuous line, and the dashed 


Mes GEpicts the run for the equilibrium initialized cases. 


PAELE TL 


Initial equilibriva aerosol values 


radius, un 


Oc Dac 5.10 10 
W 63.6 Pd Bes Bd 34:0 Oss-3 
a 9.0 10 10 520 120 
n 11 18 25 13 6.4 
d 13 20 30 15 15 
15 33 Ze 22 22 
r/s 18 28 35 28 28 
TABLE I1 
Initial observed aerosol values 
radius, um 
Oi 20 BU 10 
Betay run 42.2 Zl Zane 0 
6 May run 45.4 2230 ey eue Brest 
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Mmetecensidering the equilibriun initial v O 
Zable one, it is easy tc see the relaticnship between wind 
Seeja and froduction. In almost all cases, the aarcscl 
volume ccun*s increase aS the wind speed incr2asées. The 
distributicn of radii for a given wind speed, however, is 
feemeasc well defined. ECE BS Cwore wi ndmcpeed-, micmhe= Col 
radii are few ir nunker, and as the wind speed increases 
there is a dramatic increase in the number of larger p 
Seles. This is reasonatle due te white cap productic 
immerse appears to be a slight maximum of particles at Z.0 
mMicrens, regquardless of wind speed. The 15 micron initiali- 
Zaticn values inthe tables and the model cutput fcr the 
same size are not shown. This is because the ckserved 


meee ael values collected frem the aircraft were zero. 


ie CCtparing the MAGAT initial values, one finds that the 
smaller radii are approximately twice the equilibrium count 
for é€ven the largest wind sreed. ests tT Hteresta na ste, Nees 
that for such a larce variation in initial wind e¢peed for 
fmmer two cays, there is very little variaticn in the initial 
values, for all the the particle sizes. Based on results of 
Monahan [et al 1983/, Lewes = obcorsuaprising te find no 
volume ccunts for the larger particle sizes since wind speed 
was akove 10 m/sec. For the strong wind case, there were no 


memoscclisc abcve 10 microns collected from the aircraft. 


D. FIRST MCDEL HUN FESULTS 


mine Glcts of dv/dr over the 24 hour period for the first 
model run are shewn in figures 4.5 and 4.6. Fer each radius, 
the sclid line indicates the profile with the MAGAT initial- 
ized values, and the dashed line represents the profile with 


the equilibrium initialization. 
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Home Lhe three saaller sizes, the model shows little 
change cver the pericd, however the nodel does produce ea 
trend in aerosol proéuction. whatever volumes were use2 to 
Picaaiize, the number tends to stay within one crder of 
Megnatude cver the period. The verrttcaticon valusce c e 
very weil to the model trun with the MAGAT initializeticn. 
The medel and verification dv/dr values differ by les 
I@egnicte. For the larger size, 10 microns, note the quick 
Beoducticn of aéroscl when initialized with with 2 zero 
meesue. Within fcur tours, the model is producing cicse to 
the same values as the equilibrium case. From this feur-hcur 
point, the two cases btehave very much the same. There is a 
diurnal variation evident in the larger size cases; 4 lag is 
Beetceable in the preduction through the night time, wit! 
increase after sunrise and threugh the morning. The verifi- 
caticn values for these larger sizes are not gocd for either 


Veriafticaticn scheme. 
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For RADIUS= .8 um 
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Figures 4.5 Aerosol rlots for first run, with a) r= 0.8 ua, 
and b) r= 2.0 unm. 
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For RADIUS= 5 um 
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Figure 4.6 Same as figure 4.5, except a) r= 5.0 unm, 
and b) r= 10 um. 
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Eo SECCND MODEL RUN EESULTS 


Aerosci flots for toth the equilibrium and MAGAT initiali- 
Bea-icn schemes are shown in figure 4.7 and 4.8. Ths dv/dr 
values for both schemes are higher, due mainly to ths 
increased wind speed. For all of the particle sizes, there 
memmecre cf a diutnal variation. This is slightly e¢vidert in 
the smaller radii, and well defined in the larger radii. 
This variation is strengthened by a wind speed minimum just 
kefore sunrise. Noite Sele cyenatwene diep 14 puecdver cn, of 
mmeecne credit areund C000 PDT coincides with both the siigh= 
Gop in wind speed frcm 10 te & m/sec, and the formaticr of 
clouds, evident in the LCL and inversion figures. This 
feezease in producticn due te cloud formation is not notice- 
eee in the first model run. The two initializaticn schenes 
are ir more agreement than with the first run, again due 
Mainly tc the higher wind. With the larger wind speeds in 
this run, the model is generating vaiues closer to the veri- 
fication values (in Ecth initialization schemes) than in the 


Mees run, where the wind speed is Significantly less. 
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For RARDIUS= .8 um 
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Figure 4.7 Aercsol flots for second run, a) r= 0.8 um, and 
b) r= 2.0 un.. 
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For RADIUS= 5 um 
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Figure 4.8 Same as figure 4.7, except a) r= 5.0 um and 
b) r= 10 un. 
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Eee SUUMAR® OF RESULTS 


The mcdel is prodtcing, with a minimum of adjustment, an 
Seeurate descrigrcicn of the atmospheric boundary Ileyer 
Beeametsts. Once a suksidence rate is selected such that +t 
Beeimcur predicted irversion height approaches the cbheerved 
valué, the values cf temperature and huni 
toamemum cf adjustment. In bcth cases the verification va 
cf temperature were within 2° C, and the values of specific 
humidity within 1 gmykg. MOEeG4slOniatcart ly? the trends 
ever the 24-hour periods fer both temperature and hu 
were accurately described by the model, as déefin 


Meme ricaticn data. 


Ween re stect to the -eaeroscl =nout, the model is alse 
producing the correct trend in behavior. There is év 
mmene Cucput fcr the trend in aerosol production to be a 

oaks 


aerso the concentraticn of the initial aerosol values. This 


function of not cnly wind epeed and relative humidiz 


meme s€en in the first run, Perticie “Sezes (028 vend 2.0 
ticrcrs. Tke atmospheric parameters are the same in beth 
eaces, fcr ¢ach radits; but it is the initialization scheme 
that generates a different trend in each case. The mcdel run 
with the MAGAT initial values produces a loss in producticn, 
and the verification values are within one crder of magni- 
tude cf the model values. With the equilibrium case, the 
trend is an increasé€ in production over the time pericd, 
with the mcdel verification values two orders of magnitude 


away from the observed verification values. 


The production trends in the second mod¢el run, with both 


a 
- 
” 
Q 
m 
(n 
(D 


initialization schemes, are much the same. In 
there is less of a difference Letween the model verification 
values fer the two initializations. With the excepticn of 


the smallest radius, 0.8 microns, the verification valves 





Bemedlt the radii of the seccnd run are within one crd¢r c= 
magnitude cf the mecdel output. Tie fact “nacethe AAGAT 
3nitial values are closer tcan equilibrium v 
higher wind speed seems tec imply a less sensitive ne 
when initialized with a higher wind speed. NGCen d=. Tt 
eemeer=cl prcduction is high;. and the model accurately 
Peeaicts the behavicr of the aerosol in the stronger wind 


Case. 


oN 
(ad 





Ve CONCIUSIONS 


mee ECundary layer model, with @ minimum cf input, will 
Peoauce an accurate description of the MABL [Davidson ee ail. 
1984] and the aeroscl behavicr within the marine boundary 
Heyer. Ihe model is generating values for marine perticles 


Within this layer, and whenever there is a large diffezence 
a 


ketween the predicted and cbserved values , i. Was ene 
observed values that were toc high. Pius = Onset vapor 
suggests t S rece di. 


at the observed values collected by th a 
£ marine and ccntinental particles) may have 


h 
h 
(madé€é ufo 
included mo 


Gr® cOntinesgwal then Origianelly expected. A 
dinear distribution (vertically, through the boundary layer) 


was assumed in the aerosol initialization process. 


Mcre significantly, the model is correctly identifying 


the production trend cf aerescls over the time period. The 
model is run with identical atmospheric inputs, and 
different aerosol initialization schemes produce two 
different production trends. The wind speed and relative 


humidity (the cnly atmospheric inputs of previous models) 
are the same, and the model is generating different prefiles 
Fesed On initial aercsecl concentrations. With a minimum cf 
two case studies, it is clear that the process of entrain- 
ment cf air from above the mixed layer, containing ne marine 
Farticles, is affecting the behavior of the model output. 
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